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Abstract 

Linear covariant gauges, such as Feynman gauge, are very useful in perturbative 
calculations. Their nonperturbative formulation is, however, highly non-trivial. In 
particular, it is a challenge to define linear covariant gauges on a lattice. We con- 
sider a class of gauges in lattice gauge theory that coincides with the perturbative 
definition of linear covariant gauges in the formal continuum limit. The correspond- 
ing gauge-fixing procedure is described and analyzed in detail, with an application 
to the pure SU(2) case. In addition, results for the gluon propagator in the two- 
dimensional case are given. 
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1 Introduction 



Many direct experimental tests of QCD are performed at large energies. In this regime, 
the contributions obtainable in perturbation theory are dominant. Such perturba- 
tive calculations are often most conveniently done in linear covariant gauges, such 
as Feynman gauge (see for example [lj). On the other hand, the continuation of these 
gauges to the non-perturbative regime and, in particular, their definition on the lat- 
tice is not trivial. Only for one particular representative of this class of gauges, i.e. 
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the Landau gauge, have non-perturbative studies on the lattice been thoroughly car- 
ried out. For this gauge, progress has been made in the understanding of the infrared 
sector of the theory, in particular through study of correlation functions. This was 
possible by a combination of various methods, such as Dyson- Schwinger equations 
[2H3R .5 6 7 8 9]nUlllllll2|ll3|ll4|ll5j . renormalization-group methods [TTfl6|ll7j . stochastic 
quantization [T81IT9] . lattice calculations [20lf2Tlf22lf23lf2^U25lf26lf27Tf2^ 
and studies based on the Gribov-Zwanziger action [15]. 

In covariant gauges beyond Landau gauge, only a few results are presently available 
[46,47,48,49,50). Some of these [48,49j suggest that the Gribov-Zwanziger confinement 
mechanism [T8lll9f51f52ll53] . proposed for Coulomb and Landau gauge, may also be at 
work in the complete class of linear covariant gauges. Nevertheless, a full understanding 
of how non-perturbative effects, and especially confinement, manifest themselves in 
linear covariant gauges is still lacking. 

The aim of this work is to provide a non-perturbative implementation of linear covariant 
gauges in lattice gauge theory that goes beyond previous attempts [IB] . In Section [2] 
the non-perturbative definition of linear covariant gauges will be discussed. In Section 
[3] one possible implementation on the lattice will be introduced. The properties of the 
proposed algorithm will be studied in Section HJ showing that this gauge definition 
tends to the definition of the ordinary linear covariant gauge in the continuum limit. 
Finally, as an example, the gluon propagator is presented in Sections 14.21 and [51 We 
consider the two and three-dimensional cases and the SU(2) gauge group in order 
to study quantitatively the approach to the continuum. Our results are summarized 
in Section [6j Technical details for the generation of configurations using the Wilson 
action and for the Landau gauge fixing can be obtained from Ref. [30J. The sets of 
configurations (in the 2d case) employed for this work are listed separately for each of 
the three considered physical volumes in Tables [T][31 



2 Non-perturbative definition of linear covariant gauges 

Linear covariant gauges are defined by the average over gauge configurations satisfying 
d.A^x) = A a (x) , (1) 



for arbitrary real-valued functions A a (x). [Here A a ^(x) is the gluon field and a is the 
color index, taking values 1,2,..., iV| — 1 for the SU(iV c ) gauge group.] This average 
is usually done by considering the partition function 



Z = J PAexp j-1 J d d x £ [A a (x)f j Z(A) 



(2) 
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V 1 / 2 [fm] 


N 


P 


1/a [GeV] 


a [fm] 


Sweeps 


Configs 


3.55 


20 


10 


1.11 


0.178 


30 


2623 


3.55 


30 


21.95 


1.66 


0.118 


40 


2534 


3.55 


40 


38.6 


2.22 


0.0888 


50 


2927 


3.55 


50 


60 


2.77 


0.0710 


60 


4654 


3.55 


60 


86.3 


3.32 


0.0592 


70 


3563 


3.55 


70 


117.4 


3.88 


0.0507 


80 


4275 


3.55 


80 


153 


4.43 


0.0444 


90 


6107 


3.55 


90 


193.7 


4.99 


0.0394 


100 


8328 


3.55 


100 


239 


5.54 


0.0355 


110 


7195 


3.55 


120 


344 


6.65 


0.0296 


130 


7580 


3.55 


140 


468 


7.76 


0.0254 


150 


5227 


3.55 


160 


611 


8.86 


0.0222 


170 


2982 


3.55 


180 


773 


9.97 


0.0197 


190 


2204 


3.55 


200 


954 


11.1 


0.0178 


210 


1032 


3.55 


220 


1155 


12.2 


0.0162 


230 


959 


3.55 


242 


1397 


13.4 


0.0147 


252 


587 


3.55 


266 


1689 


14.8 


0.0133 


276 


210 


3.55 


294 


2062 


16.3 


0.0121 


304 


166 


3.55 


324 


2505 


18.0 


0.0110 


334 


87 



Table 1 

Sets of configurations used for this work for the smallest physical volume in the 2d case. 
The lattice spacing a has been evaluated using the exact infinite-volume result for the string 
tension [5l] and the input value yfa = 440 MeV. N is the lattice size in lattice units, i.e. 
iV = V 1 / 2 /a. Sweeps indicates the number of hybrid-overrelaxation sweeps between two 
consecutive thermalized configurations. Configs represents the number of thermalized config- 
urations generated in each case. 



i.e. by using a Gaussian average of width y/£ with the gauge-fixing parameter £. (Note 
that, with the usual convention h — c — 1, this gauge- fixing parameter is dimensionless 
for any space-time dimension d.) At the same time, the mean of the distribution is the 
(in perturbation theory unique) Landau-gauge configuration 

d,A*(x) = 0. (3) 
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V 1 / 2 [fm] 


N 


P 


1/a [GeV] 


a [fm] 


Sweeps 


Configs 


7.11 


40 


10 


1.11 


0.178 


50 


2633 


7.11 


50 


15.35 


1.39 


0.142 


60 


2916 


7.11 


60 


21.9 


1.66 


0.119 


70 


3973 


7.11 


70 


29.63 


1.94 


0.102 


80 


4263 


7.11 


80 


38.55 


2.22 


0.0889 


90 


4466 


7.11 


90 


48.65 


2.49 


0.0790 


100 


6507 


7.11 


100 


59.9 


2.77 


0.0711 


110 


6534 


7.11 


120 


86.1 


3.32 


0.0593 


130 


6785 


7.11 


140 


117 


3.88 


0.0508 


150 


5874 


7.11 


160 


152.7 


4.43 


0.0444 


170 


3950 


7.11 


180 


193 


4.99 


0.0395 


190 


1744 


7.11 


200 


238.4 


5.54 


0.0356 


210 


1391 


7.11 


220 


288 


6.09 


0.0323 


230 


605 


7.11 


242 


349 


6.72 


0.0294 


252 


457 


7.11 


266 


421 


7.38 


0.0267 


276 


374 


7.11 


294 


514.6 


8.09 


0.0242 


304 


134 


7.11 


324 


625 


8.99 


0.0219 


334 


80 



Table 2 

Same as in Table [H but for the intermediate physical volume. 

Thus, in perturbation theory, the linear covariant gauge is an average over the complete 
gauge orbit [as A a (x) can take any value], with Gaussian weight around the Landau 
gauge. 

The problem to be solved is to introduce a non-perturbative gauge prescription and 
to discretize it on a lattice. Of course, at the non-perturbative level, effects due to 
Gribov copies may be difficult to single out and to quantify. Nevertheless, in the lattice 
gauge-fixing procedure introduced in the next section, the determined gauge copy is 
actually unique, up to ambiguities involved in defining the (minimal) Landau gauge 
and to (possible) numerical artifacts. 

Let us note that previous attempts to formulate this class of linear covariant gauges 
[4~6] relied on the absence of (non-trivial) zero-modes of the Faddeev-Popov operator 
on the lattice in the Landau gauge. On the other hand we know that, in the infinite- 
volume limit, Landau-gauge configurations belong to the so-called first Gribov horizon 



4 



V 1 / 2 [fm] 


N 


P 


1/a [GeV] 


a [fm] 


Sweeps 


Configs 


14.2 


80 


10 


1.11 


0.178 


90 


4111 


14.2 


90 


12.55 


1.25 


0.158 


100 


4736 


14.2 


100 


15.4 


1.39 


0.142 


110 


5455 


14.2 


120 


21.9 


1.66 


0.118 


130 


6724 


14.2 


140 


29.7 


1.94 


0.101 


150 


6245 


14.2 


160 


38.6 


2.22 


0.0888 


170 


4372 


14.2 


180 


48.8 


2.49 


0.0789 


190 


2108 


14.2 


200 


60.1 


2.77 


0.0710 


210 


1533 


14.2 


220 


72.6 


3.05 


0.0647 


230 


582 


14.2 


242 


87.8 


3.35 


0.0588 


252 


596 


14.2 


266 


103 


3.69 


0.0534 


276 


465 


14.2 


294 


129.4 


4.08 


0.0483 


304 


119 


14.2 


324 


157 


4.50 


0.0438 


334 


107 



Table 3 

Same as in Table [H but for the largest physical volume. 

[30|34|55|f56] . i.e. such zero-modes do exist. Thus, it is not clear whether the procedure 
described in [46J is well defined. Moreover, in that case one obtains on the lattice a 
Faddeev-Popov operator which is different from the corresponding continuum operator. 
Here we use a different approach, which is not based on the above assumption. This 
gauge-fixing procedure is described in the following section and it will be called a 
quasi-linear covariant gauge. 



3 Quasi-linear covariant gauges on the lattice 

Our gauge-fixing prescription is a very direct implementation of the prescription for 
performing a Gaussian average around the Landau gauge. 

The first step is to transform a given (thermalized) gauge-field configuration to Landau 
gauge (see for example [oT|58|59||dD] ). Note that, in principle, in this step it is necessary 
to fix the Landau gauge completely, i.e. to select a configuration within the fundamental 
modular region [52ll6Tj . corresponding to the absolute minimum of the minimizing 
functional. However, the practical implementation of this step is still an open problem 
[27f37] . On the other hand, it has been conjectured [62] that, if one determines only 
correlation functions of a finite number of field operators, it is sufficient to fix the gauge 



5 



to a copy inside the first Gribov region, i.e. these Gribov copies should be equivalent 
for very large volumes. Hence, in our quasi-linear covariant gauge, the Landau gauge 
configuration is chosen inside the first Gribov region. 

The second step is to obtain a configuration that satisfies Eq. (pQ), for a given A a (x), 
starting from the configuration fixed to Landau gauge (in the first step). This can be 
obtained by recalling that, for an infinitesimal gauge transformation 

g{x)^l + i^\x)r\ (4) 

the gluon field A a ^(x) gets modified to 

+ (5) 



Here we indicate with 1 the identity matrix, r b are the generators of the SU(iV c ) Lie 
algebra, D ab is the covariant derivative defined as 

Df^d, + g f abc A^x), (6) 



g is the bare coupling constant and f abc are the structure constants of the gauge group. 
Thus, if one finds solution of the equations 

0„ (d^) (x) = A a (x) , (7) 



then we have 



0, 



(a; + d 



ab 



(X 



A a (x) , 



(8) 



where we used the fact that the original gluon field A a ^(x) satisfies the Landau-gauge 
condition d^A^x) = 0. Of course, this procedure is correct only if the gauge transfor- 
mation is small. On the other hand, as we will see below, terms that are of higher order 
in <ft b (x) should become smaller in the continuum limit. Let us also note that Eq. ([7]) 
above can easily be solved, using for example a Conjugate Gradient iterative method, 
since the Landau-gauge Faddeev-Popov matrix M = —d^D ab is semi-positive definite. 
Of course, for a finite lattice volume V, the functions A a (x) should be orthogonal to 
the trivial (constant) zero-modes of —d^D ab . Numerically this means that the constant 
mode has to be explicitly removed from A a (x) before starting the inversion of M. In the 
infinite- volume limit the functions A a (x) should also be orthogonal to the non-trivial 
zero-modes3 of M. 

1 We did not investigate if the presence of these non-trivial zero- modes hampers the definition 
of linear covariant gauges for £ ^ 0. 
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Furthermore, on the lattice, the gauge transformation is applied to the group-valued 
link variables U^(x) and for a given a (x) we need to define a group- valued gauge 
transformation g(x). In the SU(2) case this can easily be done by considering 



. , id) sm „ 

g = lcos(cj>)+ , <r a (9) 



= V^°, (10) 

where a a are the Pauli matrices and summation over color indices is understood!!] 
Clearly, in the limit of small <fi one finds the infinitesimal gauge transformation (jlj) , if 
we identify the Pauli matrices with the generators r a of the SU(2) Lie algebra. For finite 
a (x), the algebra- valued fields A' a (x), obtained from the link variables after this gauge 
transformation, will thus not satisfy exactly the gauge condition (jSJ). Nevertheless the 
agreement is correct up to corrections of order 0(a), becoming exact in the formal 
continuum limit. Indeed, besides the usual 0(a 2 ) errors, induced by the definition of 
the gluon field, Eq. (jjj) is correct only for the infinitesimal gauge transformation (j4]). By 
considering the finite gauge transformation ([9]), one indroduces errors of order 0(<f) 2 ). 
From Eqs. (jSJ) and (jSJ) one sees that <fi has dimension 1/go, i.e. it has mass dimension 
d/2 — 2, where d is the space-time dimension. Note that, to be more consistent, one 
should write the infinitesimal gauge transformation (jSJ) as g(x) ~ l + igodj (x)r . Thus, 
in the continuum limit these errors go indeed to zero. Of course, it would be interesting 
to study possible improvements for Eq. ([7]), in order to speed-up the approach to the 
continuum. This also shows that our approach is not unique and other choices for the 
gauge transformation g(x) are possible, leading to different sub leading behaviors. 

Finally, in the third step, we need to average over gauge copies satisfying (JTJ). As said 
above, the quantities A a (x) are randomly chosen using a Gaussian distribution with 
null mean value and width a/£. This implies that A a (x) is unbounded. On the other 
hand, on the lattice, using a compact formulation, the field A a ^(x) is bounded and so 
is its derivative d^A^(x). Thus, d^A®(x) cannot really obey a Gaussian distribution 
and in general it is not possible to satisfy Eq. JTJ on a discrete lattice for an arbitrary 
function A a (x) (see Appendix B.2 in [S3])- Nevertheless, we can still define a gauge- 
fixing procedure that becomes a linear covariant gauge in the formal continuum limit. 
To this end we note that, on the lattice, the Gaussian distribution in Eq. (J2J) can be 
written as 



J DAexp 



2£ 7i I J 2« 



2 This procedure can in principle be extended to any simple Lie group. In particular, the 
determination of the generators <p a of the gauge transformation g, using the prescription (JTJ), 
is valid for any gauge group. On the other hand, a simple map relating <p a to g, such as Eq. 
©, is usually not available and one should probably rely on numerical methods in order to 
obtain the group element g. 
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n-dependence of the transverse gluon propagator 



| n-dependence of the longitudinal gluon dressing function | 
0.07 




2.5 3 

p [GeV] 



2.5 3 

p [GeV] 



Fig. 1. The transverse gluon propagator (left panel) and the longitudinal gluon dressing 
function (right panel) as a function of p for several values of the number of samples n on 
each gauge orbit. Open circles correspond to n = 2, closed circles to n = 16, open squares 
to n = 32 and open stars to n = 1024. Results are for a 40 2 lattice at f3 = 10 for 678, 
528, 480, and 228 configurations, respectively. Note that the error bars are determined by 
considering each representative on a gauge orbit as an individual measurement. For details 
on this evaluation of the gluon propagator see Sections 14.21 and 



Here, the factor (3/(2N c ) is necessary in order to obtain the correct continuum limit. 
Indeed, working in the continuum in a generic <i-dimensional space, Eq. ([1]) can be 
made dimensionless by multiplying both sides by a 2 go. By recalling that in d dimensions 
(3 = 2iV c /(a 4 ~ d go) ^ * s clear that the lattice quantity 



becomes 
1 1 



d d x 



E Km 6 



. x 



1 

2? 



d d x 



x 



(13) 



in the formal continuum limit. Thus, on the lattice, the function A b (x) is generated from 
a Gaussian distribution with width a = J 2N c t; / [3 at each site and for each color. Note 

that in the continuum limit j3 — > oo the width of the distribution J 2N C £ / (3 shrinks to 
zero. Therefore, the functions A b (x) depart less from zero and the gauge fixing can be 
achieved using a small gauge transformation for nearly all gauge copies entering into 
the average ([2]). Note that this process is slower for larger £, since the width of the 
Gaussian scales with 



S 



Finally, if one considers the partition function ([2]) it is clear that, in order to evaluate 
an observable in the linear covariant gauge, one needs to perform two averages: one is 
the usual average over configurations C, while the other is the one along the gauge orbit 
Q. Thus, the space to be sampled is C x Q and errors should be estimated by regarding 
each sample from C x Q as an independent configuration. Of course, in order to sample 
correctly this product space, one needs to generate enough thermalized configurations 
and, for each configuration, a sufficient number n of Gaussian-distributed copies on the 
given gauge orbit. How large n should be clearly depends on the considered observables. 
For correlation functions such as the gluon propagator we find that, if the number of 
thermalized configurations is large enough, then a very small number for n yields a 
result that does not change significantly by increasing n, as shown in Fig. [TJ This 
implies that the fluctuations along a gauge orbit are smaller than between different 
orbits. 



4 Properties of the gauge-fixing 

In this section we study the extrapolation to the continuum limit of the gauge-fixing 
condition described in the previous section. The results depend, of course, on the value 
of £, i.e. on the width of the Gaussian distribution. Here we consider the case £ = 1 
(Feynman gauge) for the calculations in three dimensions. Below, in two dimensions, 
we will consider the case £ = 1/100. 

4-1 First moment 

Our first check refers to the distributions of the one-point correlation function A' a (:r) = 
dpA! a (x), where A' a {x) is defined in Eq. (j5J). Since A a (x) is generated using a Gaussian 
distribution, one should in principle obtain that the quantity A /a (a;), evaluated after 
fixing the configuration to the quasi-linear covariant gauge, should also obey a Gaussian 
distribution with the same width. On the other hand, due to discretization errors and 
to the finite-precision arithmetic employed in the determination of <p a (x), one can 
expect deviations from the original Gaussian distribution. Thus, this check is a simple 
and direct way to measure the quality of the gauge fixing and the deviation from the 
continuum gauge condition (jSJ). Results are reported in Fig. [2] and in Table 0] for three 
different lattice volumes in 3d with (roughly) the same physical volume (2.1 fm) 3 . We 
find that in the three cases the width of the distribution of A' a (x) is in rather good 
agreement with the original distribution, since the measured value £ m of the gauge- 
fixing parameter is very close to 1. The deviation from 1 is actually expected, and 
corresponds to the necessity of renormalizing £ in linear covariant gauges. We also 
find that the gauge-fixing procedure does not introduce skewness, while in the three 
cases a non-vanishing kurtosis is clearly observable. On the other hand, this kurtosis is 
decreasing with increasing j3 (see Table 0J. Also note (see Fig. [2]) that the distribution 
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Distribution at (3=4.2 



io*r 
10 - 




Distribution at (3=6.0 



10 7 
10 6 
10= 
10 4 
10 3 
10 2 
10 

1 





Fig. 2. The measured distribution of A' a (a;) (solid line), for all lattice sites x, for (3 = 4.2 and 
V = 12 3 (top left panel), = 6.0 and V = 18 3 (top right) and = 9.0 and V = 28 3 (bottom). 
In the three cases the physical volume is roughly (2.1 fm) 3 . In all cases we used 10 different 
(thermalized) configurations with n = 256 copies on each gauge orbit (we performed 200 
thermalization sweeps for the initial configurations and 29 sweeps between two configurations, 
using the standard hybrid-overrelaxation-scheme [30]). The dashed line represents the initial 
(Gaussian) distribution with width y/2N c £/(3. The plot range along the x-axis indicates the 
range of values obtained for A' a (x). 



p 


N 3 




Cm 


Kurtosis 


4.2 


12 3 


1 


0.978 


0.539 


6.0 


18 3 


1 


0.953 


0.124 


9.0 


28 3 


1 


0.951 


0.00977 



Table 4 

Properties of the measured distribution of A' a (x) 



of A /a (x) becomes visibly different from that of A a (x) for values on the y-axis of about 
10 4 at = 4.2, 10 3 &tp = 6.0 and 10 2 at p = 9.0. 
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Relative deviation at p=4.2 | Relative deviation at p=6.0 | 




Fig. 3. Histogram of the relative pointwise deviation 5 for the same sets of configurations 
considered in Fig. [2j Results are shown for (3 = 4.2 in the top left panel, for (3 = 6.0 in the 
top right panel and for (3 = 9.0 in the bottom left panel. The x-axis has been truncated at 
5 = 10 4 . For comparison, the situation in Landau gauge according to (115p is shown. The 
result is for a 12 3 lattice at f3 = 4.2, and 5 has been truncated at 10~ 13 with the largest 5 at 
around 10~ 12 . All distributions have been normalized to a unit area. 



We also evaluated the relative (pointwise) deviation 

(14) 



|A a (x) - cVWI 



|A-(x)| 



Results are shown in Fig. [3] for the same set of configurations considered in Fig. [2J 
Clearly we always find sites where S is very large. However, the maximal value of 5 as 
well as the number of sites with large relative deviation, both in absolute number and 
normalized to the total number of sites, is decreasing with increasing (3, even though 
this is not clearly visible from the plots. E. g., at the cutoff value S = 10 4 considered in 
Fig. [31 the value in the y-axis is about 6 x 10~ 6 , 4 x 10~ 6 and 3 x 10~ 6 , respectively at 
(3 = 4.2, 6.0 and 9.0. Since the histograms are normalized to a unit area, this implies 
that the probability of such a large deviation is reduced at large (3. A possible source 
for the very large deviations is likely related to functions A a that exceed the range that 
can be covered by and d^A^ using the compact formulation on the lattice. 
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We have also checked that in Landau gauge the distribution for 

6=\d,A;(x)\ 



(15) 



is quite similar to those shown in Fig. [3j although the absolute size of 5 is then much 
smaller and of the order of the quality of the gauge-fixing. This is also displayed in 
figure [31 



4-2 Second moment: The longitudinal gluon propagator 



We have seen in the above section that, if one considers single-site properties, the 
quasi-linear covariant gauge indeed seems to approach the usual linear covariant gauge 
in the limit (3 — > oo. On the other hand, we should check if the same is true for 
all correlation functions of d^A^x). Here we consider only the two-point correlation 
functions of d^A a (x), since higher correlation functions are statistically very noisy |30j . 

In the case of the two-point correlation function 

<^{x)dlA h M> (16) 



the distribution should be proportional to 8 ab 8 tiv 8(x — y). By applying the Fourier- 
transform one finds the well-known Slavnov- Taylor identity (STI) that the longitudinal 
part of the gluon propagator does not become dressed and is proportional to 5 ab /p 2 , 
with the renormalized gauge parameter £ r appearing in the constant of proportionality 

DP 

p,p u Df v {p) = i r 5 ab . (17) 



Let us recall that this STI is a trivial consequence of the gauge condition. Thus, its 
violation as a function of (3 is a useful check in order to control the extrapolation to 
the continuum limit. 

Results for the longitudinal part of the gluon propagator, defined a 



DM ^ kwtw^ < 18) 



3 Here we did not consider the mid-link Fourier transform of the gluon field [64J. However, 
the behaviors of the longitudinal gluon propagator for momenta along the x-axis and along 
the xt-diagonal are rather similar (see Fig. H|), suggesting that discretization effects related 
to the definition of the gluon field in momentum space are quite small. 
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Longitudinal qluon propagator [(14.2 fm) 2 1 If Longitudinal gluon dressing function [(14.2 fm) 2 ] \ 




4.5 5 
p [GeV] 



Longitudinal gluon propagator [(3.56 fm) 2 1 II Longitudinal gluon dressing function [(3.56 fm) 2 ]~| 



30 : 



10- 



10" 




o 



3Jb 



I I I I I I I I I I I I I I I I I I I I I I I I I I M I I I I I I I I 



0.5 1 1.5 2 2.5 3 3.5 



p [GeV] 



4 4.5 5 
p [GeV] 



Fig. 4. The longitudinal part of the gluon propagator for three different lattice volumes V: 
(14.2 fm) 2 (top), (7.11 fm) 2 (middle) and (3.56 fm) 2 (bottom). Symbols refer to different (3 
values in the three different panels, i.e. (3 = 48.8 (Circles), 87.8 (Squares), 129.4 (Triangles) in 
the top panel; (3 = 59.9 (Circles), 117 (Squares), 238.4 (Triangles) in the middle panel; (3 = 10 
(Circles), 1155 (Squares), 2505 (Triangles) in the bottom panel. All dressing functions have 
been renormalized such that their lowest value coincides with the renormalized continuum 
value £ = 1/100, indicated by a line. In all cases we considered n=16 gauge copies for each 
gauge orbit. Momenta are along the x-axis and along the xt-diagonal. 
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are reported in Fig. |4]for the case £ = 1/100. This value of £ was chosen, instead of the 
Feynman gauge £ = 1, in order to speed up the approach to the continuum limit. One 
clearly sees that the dressing function p 2 D(p) is far from flat, but is becoming slowly 
flatter with increasing (3. Only at large momenta does one find an essentially constant 
function. However, also at large momenta, renormalization effects are clearly visible. 
They can be quantified by evaluating the renormalization factors using the relation 
and renormalization condition 

^ r = Z i U = i- (19) 



Results are reported in Table [5J The increase of with (3 is consistent with a log- 
arithmic evolution, in accordance with the perturbative expectation. Moreover, this 
evolution depends on the considered physical lattice volume, being faster for larger 
volumes. As Z^ tends to infinity in the limit a — > 0, the attraction of the Landau gauge 
as a fixed point is evident. 

In any case, the main finite-/? effects can be seen at low momenta. Indeed, there is 
a clear maximum in the longitudinal gluon dressing function for momenta of about 
0.5-1 GeV (see Fig. H]). We also find that this maximum seems to become flatter as 
(3 increases, even though this appears as a remarkably slow process. In addition, the 
height of the maximum increases with the physical volume. The j3 dependence of this 
maximum is shown in Fig. [51 One sees that, for small (3, the maximum seems to increase 
first, before starting to decrease slowly to the expected continuum value PI The slow 
evolution with f3 is probably related to the fact that (3 enters only as a square-root in the 
width a; it may also be influenced by the logarithmic running of renormalization effects 
in two dimensions. The increase in strength of this artifact as the volume increases is 
possibly due to the development of would-be-zero modes of the Faddeev-Popov operator 
with volume. Projecting out these modes from A a could reduce the strength of these 
discretization artifacts. 

It is also interesting to observe from Fig. [5] that the points characterized by very large 
errors are also those that deviate most from the expected behavior. The reason for 
this is a strongly asymmetric distribution of the longitudinal part of the gluon dressing 
function as a function of configuration. Particularly large values may be due to ex- 
ceptional configurations [30J. At the same time, this could explain the large statistics 
required to obtain an essentially continuous function of j3 in Fig. [51 Also, due to com- 
putational limitations, we did not check if these results can be improved by increasing 
n at large (3. 

Finally, we find similar qualitative results when considering a space-time dimension d 
larger than two. 



4 This kind of non-monotonic behavior has also been observed in the 2d Landau case when 
considering finite- volume effects |34| . 



14 



N 


P 


L [fin] 




P 


L [fin] 




P 


L [fm] 




20 


10 


3.55 


1.003(2) 














30 


21.95 


3.55 


1.010(2) 














40 


38.6 


3.55 


1.013(2) 


10 


7.11 


1.030(7) 








50 


60 


3.55 


1.013(1) 


15.35 


7.11 


1.00(1) 








60 


86.3 


3.55 


1.014(3) 


21.9 


7.11 


1.031(3) 








70 


117.4 


3.55 


1.012(2) 


29.63 


7.11 


1.045(3) 








80 


153 


3.55 


1.013(1) 


38.55 


7.11 


1.053(2) 


10 


14.2 


1.017(8) 


90 


193.7 


3.55 


1.011(3) 


48.65 


7.11 


1.061(3) 


12.55 


14.2 


1.02(1) 


100 


239 


3.55 


1.011(1) 


59.9 


7.11 


1.064(1) 


15.4 


14.2 


1.088(8) 


120 


344 


3.55 


1.011(1) 


86.1 


7.11 


1.062(3) 


21.9 


14.2 


1.203(4) 


140 


468 


3.55 


1.014(1) 


117 


7.11 


1.064(1) 


29.7 


14.2 


1.231(5) 


160 


611 


3.55 


1.018(2) 


152.7 


7.11 


1.067(2) 


38.6 


14.2 


1.248(3) 


180 


773 


3.55 


1.018(2) 


193 


7.11 


1.073(3) 


48.8 


14.2 


1.21(1) 


200 


954 


3.55 


1.024(3) 


238.4 


7.11 


1.072(3) 


60.1 


14.2 


1.272(8) 


220 


1155 


3.55 


1.022(3) 


288 


7.11 


1.079(4) 


72.6 


14.2 


1.277(6) 


242 


1397 


3.55 


1.026(4) 


349 


7.11 


1.083(9) 


87.8 


14.2 


1.294(5) 


266 


1689 


3.55 


1.037(7) 


421 


7.11 


1.078(6) 


103 


14.2 


1.306(6) 


294 


2062 


3.55 


1.038(8) 


514.6 


7.11 


1.093(9) 


129.4 


14.2 


1.31(1) 


324 


2505 


3.55 


1.05(1) 


625 


7.11 


1.12(1) 


157 


14.2 


1.31(3) 



Table 5 

The renormalization constant denned in Eq. (|19|) . obtained by requiring the lowest value 
of the longitudinal dressing function to coincide with £ = 1/100. Here we did not check for 
possible Z(2) effects, which strongly affect the Landau gauge gluon propagator for large (3 
values [65] . 



5 The transverse gluon propagator 

As a first application of the gauge-fixing procedure introduced above, we consider the 
numerical determination of the transverse part of the gluon propagator. There are pre- 
dictions [48 49J that, at least in the 4d case, this propagator should be suppressed in the 
infrared limit. Results in the 2d case for various volumes and f3- values, using a standard 
definition of the transverse gluon propagator [30] . are shown in Fig. [6] for £ = 1/100, 
considering the same sets of configurations used for the longitudinal gluon propaga- 
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(3-dependence of the maximal value in the longitudinal gluon dressing function 
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Fig. 5. The /3-dependence of the maximum of p 2 Dl(p) of the (renormalized) longitudinal 
gluon propagator. Circles refer to the physical volume (3.55 fm) , squares to (7.11 fm) 2 and 
triangles to (14.2 fm) 2 (see Table [5] for the corresponding (3 values). 



tor reported in Fig. HI We see that for momenta above about half a GeV there is no 
pronounced volume or /3-dependence, in contrast to the case of the longitudinal propa- 
gator. In particular, the dressing function is qualitatively similar to the corresponding 
dressing function in Landau gauge [31]. On the other hand, at small momenta (i.e. 
below half a GeV) we find finite-volume and /3-effects. In particular, similarly to the 
Landau-gauge case, there seems to be a maximum in the gluon propagator for momenta 
of about 400 MeV, at least at large enough physical volume. However, on the largest 
physical volume, there is also a rise in the infrared (below about 250 MeV), which 
diminishes with increasing /3. Thus, the quantitative properties of this rise are strongly 
affected by discretization effects, and a final conclusion on the infrared behavior cannot 
yet be drawn. 
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Transverse gluon propagator [(14.2 fmfl | Transverse gluon dressing function [(14.2 fmf] | 
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Fig. 6. The transverse part of the gluon propagator for three different lattice volumes V: 
(14.2 fm) 2 (top), (7.11 fm) 2 (middle) and (3.56 fm) 2 (bottom). Symbols refer to different (3 
values in the three different panels, i.e. (3 = 48.8 (Circles), 87.8 (Squares), 129.4 (Triangles) 
in the top panel; (3 = 59.9 (Circles), 117 (Squares), 238.4 (Triangles) in the middle panel; j3 = 
10 (Circles), 1155 (Squares), 2505 (Triangles) in the bottom panel. Propagators and dressing 
functions have not been renormalized. In all cases £ = 1/100 and we considered n=16 gauge 
copies for each gauge orbit. Momenta are along the x-axis and along the xt-diagonal. 
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6 Summary 

We have described a new class of gauges on the lattice, which become the linear co- 
variant gauges in the formal continuum limit. The residual 0(a) errors have been 
investigated in detail. We find that our gauge approaches the correct continuum gauge, 
albeit slowly (especially at large physical volumes). One should check if it is possible 
to speed up the approach to the continuum by eliminating the leading discretization 
errors. This speed-up would probably be crucial in the study of the infrared sector of 
the theory, in particular when considering space-time dimension d larger than two. 

As a first application we have presented a preliminary study of the transverse part 
of the gluon propagator in two dimensions. Our results in the infrared limit are not 
conclusive and more investigation is required in order to obtain a clearer picture. Of 
course it would also be interesting to extend this study to higher dimensions and to 
other quantities, such as the ghost propagator. However, since the gauge field is no 
longer transverse, this study is more involved than in the usual Landau gauge. Finally, 
one should also check for possible Gribov-copy effects on the evaluated quantities. 
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